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ABSTRACT 

We present a physically motivated model for the early co-evolution of massive spheroidal 
galaxies and active nuclei at their centers. Within dark matter halos, forming at the rate pre- 
dicted by the canonical hierarchical clustering scenario, the gas evolution is controlled by gravity, 
radiative cooling, and heating by feedback from supernovae and from the growing active nucleus. 
Supernova heating is increasingly effective with decreasing binding energy in slowing down the 
star formation and in driving gas outflows. The more massive proto-galaxies virializing at earlier 
times are thus the sites of the faster star-formation. The correspondingly higher radiation drag 
fastens the angular momentum loss by the gas, resulting in a larger accretion rate onto the central 
black-hole. In turn, the kinetic energy carried by outflows driven by active nuclei can unbind 
the residual gas, thus halting both the star formation and the black-hole growth, in a time again 
shorter for larger halos. For the most massive galaxies the gas unbinding time is short enough 
for the bulk of the star-formation to be completed before type la supernovae can substantially 
increase the Fe abundance of the interstellar medium, thus accounting for the a-enhancement 
seen in the largest galaxies. The feedback from supernovae and from the active nucleus also de- 
termines the relationship between the black-hole mass and the mass, or the velocity dispersion, of 
the host galaxy, as well as the black-hole mass function. In both cases the model predictions are in 
excellent agreement with the observational data. Coupling the model with GRASIL (Silva et al. 
1998), the code computing in a self-consistent way the chemical and spectrophotometric evolution 
of galaxies over a very wide wavelength interval, we have obtained predictions in excellent agree- 
ment with observations for a number of observables which proved to be extremely challenging for 
all the current semi-analytic models, including the sub-mm counts and the corresponding rcdshift 
distributions, and the epoch-dependent K-band luminosity function of spheroidal galaxies. 

Subject headings: galaxies: elliptical and lenticular, lcD — galaxies: evolution — galaxies: formation — 
QSOs: formation 



1. Introduction 

Although the traditional approach to galaxy 
formation and evolution regards nuclear activity 
as an incidental diversion, it is becoming clear, 
beyond any reasonable doubt, that the formation 
of super-massive black holes powering nuclear ac- 
tivity is intimately linked to the formation of its 
host galaxy and plays a key role in shaping its 
evolution. Evidences supporting this view include: 
the discovery that Massive Dark Objects (MDOs), 
with masses in the range ~ 10 6 -3 x 10 9 M Q and 
a mass function matching that of baryons ac- 
creted onto black holes during the quasar activity 
(Salucci et al. 1999), are present in essentially all 
local galaxies with a substantial spheroidal com- 
ponent (Kormendy & Richstone 1995; Magorrian 
ct al. 1998; van der Marel 1999; Kormendy & 
Gebhardt 2001); the tight correlation between the 
MDO mass (Mbh) and the velocity dispersion of 
stars in the host galaxy (Magorrian et al 1998; 
Ferrarese & Merritt 2000; Gebhardt et al. 2000; 
Trcmaine et al. 2002), the mass of the spheroidal 
component (McLure & Dunlop 2002; Dunlop et 
al. 2003; Marconi & Hunt 2003), and the mass 
of the dark halo (Ferrarese 2002); the correspon- 
dence between the luminosity function of active 
star-forming galaxies at z ~ 3, the B-band lumi- 
nosity function of quasars and the mass function 
of dark halos at the same redshift (Haehnelt et 
al. 1998, Monaco et al. 2000); the similarity be- 
tween the evolutionary histories of the luminosity 
densities of galaxies and quasars (e.g. Cavaliere & 
Vittorini 1998). Recently Shields et al. (2003) 
have found that the correlation between Mbh and 
the stellar velocity dispersion is already present at 
redshift up to ^ 3. 

As discussed by Granato et al. (2001), the mu- 
tual feedback between galaxies and quasars during 
their early evolutionary stages may be the key to 
overcome some of the crises of the currently stan- 
dard scenario for galaxy evolution. For example, 
predictions of semi-analytic models (Devriendt & 
Guiderdoni 2000; Cole et al. 2000; Somerville et 
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al. 2001; Mcnci ct al. 2002) are persistently un- 
able to account for the surface density of massive 
galaxies at substantial redshift detected by (sub)- 
mm surveys with SCUBA and MAMBO (Blain ct 
al. 2002; Scott et al. 2002) and by deep K-band 
surveys (Cimatti et al. 2002b; Kashikawa et al. 
2003), unless ad- hoc adjustments are introduced. 
The difficulty stems from the fact that the stan- 
dard CDM scenario tends to imply that most of 
the star formation occurs in relatively small galax- 
ies that later merge to make bigger and bigger 
objects. On the contrary, the data indicate that 
galaxies detected by (sub)-mm surveys are mostly 
very massive, with very high star-formation rates 
(~ lO^loyr- 1 ), at z > 2 (Dunlop 2001; Ivison 
et al. 2002; Aretxaga et al. 2003; Chapman et 
al. 2003). All these data are more consistent with 
the traditional "monolithic" scenario, according to 
which elliptical galaxies formed most of their stars 
in a single burst, at relatively high redshifts, and 
underwent essentially passive evolution thereafter. 
On the other hand, the "monolithic" scheme is in- 
adequate to the extent that it cannot be fitted in 
a consistent scenario for structure formation from 
primordial density perturbations. 

Clues on the timing of evolution of both galax- 
ies and quasars are provided by chemical abun- 
dances (e.g. Friaca & Terlevich 1998). Spec- 
troscopic observations demonstrate, even for the 
highest redshift quasars, a fast metal enrichment 
of the circum-nuclear gas (e.g. Hamann & Fer- 
land 1999; Fan et al. 2000, 2001; Frcudling et al. 
2003). Statistical studies of local E/S0 galaxies, 
hosting super-massive black-holes, show that the 
most massive galaxies are also the most metal rich, 
the reddest, and, perhaps, the oldest (Forbes & 
Ponman 1999; Trager et al. 2000a, b). Such galax- 
ies show an excess a-elements/Fe ratio, compared 
to solar (a-enhancement; Trager et al. 2000a, b; 
Thomas et al. 2002), suggestive of very intense 
but short star-formation activity. Also observa- 
tions of hosts of high redshift quasars show that 
they are as massive as expected from the local 
^QSO-^host relation, and that most of their stars 
arc relatively old, although star forming regions 
are present (Kukula et al. 2001; Hutchings et al. 
2002; Ridgway et al. 2002; Stockton & Ridgway 
2001). 

The work carried out by our group in the last 
several years, aimed at constructing physically 
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grounded models for the joint formation and evo- 
lution of quasars and spheroidal galaxies in the 
framework of the standard hierarchical clustering 
scenario (Granato et al. 2001), has shown that, to 
account for the surface density of massive galax- 
ies at substantial redshifts detected by sub-mm 
SCUBA surveys (which turns out to be far in ex- 
cess of predisctions of standard semi-analytic mod- 
els), allowing for the observed relationships be- 
tween quasars and galaxies, it is necessary to as- 
sume that the formation of stars and of the central 
black hole took place on shorter timescales within 
more massive dark matter halos. In other words, 
the canonical hierarchical CDM scheme - small 
clumps collapse first - is reversed for baryon col- 
lapse and the formation of luminous objects ( Anti- 
hierarchical Baryon Collapse scenario). This be- 
havior, attributed to the feedback from super- 
nova explosions and, for the most massive galax- 
ies, from nuclear activity, may account simultane- 
ously for evolutionary properties of quasars and of 
massive spheroidal galaxies (Monaco et al. 2000; 
Granato et al. 2001), for clustering properties of 
SCUBA galaxies (Magliocchetti et al. 2002; Per- 
rotta et al. 2003), and for the metal abundances in 
spheroidal galaxies and in bulges of later Hubble 
types hosting super-massive black holes (Romano 
et al. 2002). 

Recently, Cattaneo & Bernardi (2003), using 
the early type galaxy sample in the Sloan Digital 
Sky Survey (Bernardi et al. 2003a), investigated 
the hypothesis that quasars formed together with 
the stellar population of ellipticals, finding a con- 
sistency with the observed luminosity function of 
optical quasars. 

However, there is not, as yet, a clear under- 
standing of the physical mechanisms governing the 
interactions among the active nucleus and the host 
galaxy. Silk & Rces (1998) and Fabian (1999) 
proposed that the relationship between the black- 
hole mass and the velocity dispersion (or mass) of 
the host stellar spheroid may be the effect of the 
quasar feedback. However, in its present form, 
the proposed mechanism does not imply shorter 
star-formation times for more massive galaxies and 
therefore cannot easily explain the a-enhancemcnt 
of more massive objects. Wang & Biermann 
(1998) suggested that the formation of both el- 
liptical galaxies and super-massive black-holes at 
their centers is related to the merging of two proto- 



disks. Kauffmann & Hacnhelt (2000) and Hacn- 
helt & Kauffmann (2000) have analyzed the evo- 
lution of active nuclei and of host galaxies in the 
framework of the hierarchical clustering scenario, 
using a semi-analytic approach. In their scheme, 
the merging process determines both the evolu- 
tion of galaxies and the growth of the black-holes 
at their centers. This model predicts a rapid evo- 
lution of galaxies with redshift. Volonteri et al. 
(2002, 2003) presented a model in which most of 
the mass in BHs is assembled in accretion episodes 
triggered by merging. They found that the galaxy 
merging would leave about 10% of massive black- 
holes distributed in galactic halos and a simi- 
lar fraction of binary super-massive black-holes in 
galactic centers. Di Matteo et al. (2003) consid- 
ered the gas content in galaxies, as predicted by 
cosmological hydrodynamical simulations includ- 
ing sub-grid prescriptions for gas cooling and star 
formation (but not for BH growth) . They pointed 
out that the observed Mbh-c correlation is well re- 
produced, provided that a linear relationship be- 
tween the gas and BH masses (at z > 1) is as- 
sumed. 

In this paper we carry out an investigation of 
the physical processes driving the growth of the 
central black hole and the effect of the energy re- 
leased by the active nucleus on the surrounding 
proto-galactic gas. The corresponding prescrip- 
tions for matter outflow, gas cooling and collapse, 
star formation, etc. are implemented in a model 
for the formation and evolution of galaxies, inter- 
faced with our code computing their spectral en- 
ergy distribution as a function of their age, taking 
into account the evolution of stellar populations, of 
metal abundances, of the dust content and its dis- 
tribution (GRASIL: Silva et al. 1998, Silva 1999, 
Granato et al. 2000). Our aim is to put on a more 
solid physical basis the approach by Granato et al. 
(2001), which partly relies on empirical recipes. 

The model has been tested against the observed 
correlation of BH masses with the mass and veloc- 
ity dispersion of the host galaxies, the observed K- 
band and sub- millimeter counts and the associated 
redshift distributions, along with the chemical and 
photometric properties of the local E/S0 galaxies. 
In subsequent papers the model will be used to 
predict the time-dependent luminosity function of 
quasars up to z > 6 and of passively evolving el- 
lipticals, and their clustering properties, with ref- 
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erence to the observational capabilities of space 
missions such as SIRTF, to be launched in a few 
months, JWST, of ground millimeter telescopes 
such as ALMA, and of the forthcoming deep near- 
IR surveys from the ground. 

In Sect. 2 we describe the basic physical 
processes governing the evolution of the star- 
formation rate in spheroidal galaxies, the growth 
of central black-holes and their feedback. In 
Sect. 3 we present our main results, which are 
discussed in Sect. 4. In Sect. 5 we summarize our 
main conclusions. 

We adopt a cold dark matter cosmology with 
cosmological constant, consistent with the Wilkin- 
son Microwave Anisotropy Probe (WMAP) data 
(Bennett et al. 2003), as well as with infor- 
mation from large scale structure (Spcrgel et al. 
2003): VL m = 0.29, fl b = 0.047, and fl A = 0.71, 
Ho = 72 km s _1 , (T 8 = 0.8, and an index n = 1.0 
for the power spectrum of primordial density fluc- 
tuations. 

2. Basic ingredients 
2.1. Dark matter halos 

Following Navarro et al. (1997) and Bullock 
ct al. (2001b) we identify as a virialized halo at 
redshift z a volume of the Universe of radius r v ; r 
enclosing an overdensity A v ; r (2;) which, for a flat 
cosmology, can be approximated by 



(18tt 2 + 82x - 39x 2 ) 
Q(i) : 



(1) 



where x = fl(z) — 1 and Q,(z) is the ratio of the 
mean matter density to the critical density at red- 
shift z (Bryan & Norman 1998). The halo mass is 
then given by 

M vir = — A vir (z) j0 „(z)r3 ir , (2) 

p u (z) being the mean universal matter density. 
Useful quantities are the rotational velocity of the 
DM halo at its virial radius, V^ r = GM v i r /r v j r , 
and the equilibrium gas temperature in the DM 
potential well 



kT = — M m pK 



vir "> 



(3) 



where m p is the proton mass and \im v is the mean 
molecular weight of the gas. 



The virialized halos exhibit a universal density 
profile (Navarro et al. 1997; Bullock et al. 2001b) 
well described by 



p{r) 



Ps 



cx(l + ex) 2 ' 



(4) 



where x = r/r v ; r and c is the concentration param- 
eter. Analyzing the profiles of a large sample of 
virialized halos obtained through high-resolution 
N-body simulations, Bullock et al. (2001b) found 
that, at any redshift, c is a weak function of 
the mass c oc M^ r 013 , while, at fixed mass, c oc 
(l + z)-\ 

Numerical simulations and theoretical argu- 
ments show that the dark matter assembly in ha- 
los proceeds through a first phase of fast accretion, 
followed by a slow phase (Wechsler et al. 2002; van 
den Bosch 2002). Zhao et al. (2003) showed that 
the potential well of the halos is built up during 
the fast accretion. The subsequent slow accretion 
does not change the "identity" of the halo, though 
significantly increasing its mass. Guided by these 
results, in the following we will assume that, for 
the mass and redshift ranges we are interested in, 
the virialization epoch of DM halos coincides with 
the end of the fast accretion phase and with the 
beginning of the vigorous star formation of the 
proto-spheroidal galaxies, and that these galax- 
ies keep their identity till the present time. Our 
approach implies, in the mass and redshift range 
considered, a one to one correspondence between 
halos and galaxies, consistent with the available 
data on clustering of Lyman-break (Bullock et al. 
2002) and SCUBA galaxies (Magliocchetti et al. 
2001; Perrotta et al. 2003). 

The formation rate of massive halos (M v j r > 
2.5 x 10 11 M Q ) at z > 1.5 is approximated by 
the positive part of the time derivative of the halo 
mass function n(M v ; r ,z) (Haehnelt & Rees 1993; 
Sasaki 1994; Peacock 1999). The negative part is 
small for this mass and redshift range, consistent 
with our assumption that massive halos survive 
till the present time. Following Press & Schechter 
(1974) we write: 



n(M vir , z) = Y72~ v f( u 



dhiv 



M 2 r " v 'd\nM vh 



(5) 



where p is the average comoving density of the uni- 
verse and v — [8 c {z) / (j${M v - lr )] 2 , as being the rms 
initial density fluctuations smoothed on a scale 
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containing a mass M V - 1T , and 8 C the critical over- 
density for spherical collapse. The latter quantity 
is given by S c (z) = 5o/b(z), 5q being the present 
day critical over-density (1.686 for an Einstein de 
Sitter cosmology, with only negligible dependen- 
cies on Q TO and Qa) and b(z) the linear growth 
factor, for which we adopt the approximation pro- 
posed by Li & Ostrikcr (2002). For the function 
vf(y) we adopt the expression given by Sheth & 
Tormcn (2002), which also takes into account the 
effect of ellipsoidal collapse: 



^)=A[l + M-f]( T ) 1 /2__, ( 6 ) 

where A = 0.322, p = 0.3 and a = 0.707. 

2.2. Star-formation rate 

Recent studies on the angular momentum of 
dark matter (DM) halos suggest that a significant 
fraction of the mass has a low specific angular mo- 
mentum (Bullock et al. 2001a). Also, the results 
of high resolution simulations including gas sug- 
gest that a significant fraction of the angular mo- 
mentum of the baryons is redistributed to the dark 
matter by dynamical friction (see e.g. Navarro & 
Steinmetz 2000), although the heating of the gas in 
sub-galactic halos may increase the tidal stripping 
and the angular momentum of the gas (Mailer & 
Dekel 2002). As far as the effect of angular mo- 
mentum can be neglected, as is probably the case 
for the formation of spheroidal galaxies, the col- 
lapse time of baryons within the host dark matter 
halo, tcoii, is the maximum between the free-fall 
time, 

t dyn (r) = [3 7 r/32Gp(r)] 1 / 2 , (7) 
and the cooling time 



icooi(r) 



3 p gas (r) kT 

2 iim p Cn 2 e {r)K(T) ' 



(8) 



both computed at the virial time £ v ir- In the above 
equations p is the total matter density, p gas is the 
gas density, n e is the electron density, A(T) is the 
cooling function and C = (n 2 (r)) j (n e (r)) 2 is the 
clumping factor, assumed constant. In the follow- 
ing we adopt the cooling function given by Suther- 
land & Dopita (1993), which includes the depen- 
dence on metal abundance. 

We consider a single-zone galaxy with three 
gas phases: diffuse gas in the outer regions, with 



mass -Mi n f(i), infalling on a dynamical timcscale, 
cool gas with mass M co id(i), available to form 
stars, and hot gas with mass Mh t(i), eventually 
outflowing. At the virialization we assume that 
Mnf(^vir) = fbM v i T , with fi, = 0.16, the univer- 
sal ratio of baryons to DM (Bennett et al. 2003; 
Spergel et al. 2003). The infalling gas, initially 
at the equilibrium temperature in the DM poten- 
tial well [Eq. (3)], is transferred to the cool star- 
forming phase at a rate: 



M C oid(t) = Mi n f(t)/max[t coo \(r v ir),tdyn(r v ir)] 
The star formation rate (SFR) is given by: 



(9) 



-r 

Jo 



dM cold (r, t) 



max[t coo \(r),t dyn (r)} dr 



dr, 



(10) 

where we assume that the cold gas distribution 
still follows the DM distribution. This assumption 
is admittedly quite unrealistic since cold cloudlets 
should rather fall towards the center. But again, 
a realistic modelling of the cold gas distribution 
would be too ambitious at the present stage. 

As mentioned above, according to recent stud- 
ies (Wechsler et al. 2002, van den Bosch 2002, 
Zhao et al. 2003) while the build-up of the poten- 
tial well is rather fast, the mass of the halo goes 
on increasing over a time scale of the order of the 
Hubble time. The baryons associated with this 
slow accretion occupy a large volume so that have 
a relatively low density and feel also the heating 
effect of the supernova and quasar feedback. It is 
therefore reasonable to assume that their cooling 
time is comparable to, or longer than, the Hubble 
time, and does not partake in the star formation 
process. 

Strictly speaking, in Eq. (9) the increase of the 
cooling time due to the dilution of the hot gas 
as the cold gas drops out should be taken into 
account. The effect is, however, minor because the 
fraction of gas in the cold phase is always < 50% 
and can therefore be ignored given the exploratory 
nature of the present model. 

Although the Initial Mass Function (IMF) may 
depend on the physical properties of the gas, such 
as density and metallicity (see Eisenhauer 2001 
and references therein), it is often assumed to 
be independent of time and galaxy mass. We 
will consider in the following the IMF preferred 
by Romano et al (2002) on the basis of chemi- 
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cal abundances in local ellipticals, i.e. <1>(M) oc 
M-° A for M < 1 M© and $(M) oc M -1 ' 25 for 
1M Q < M < 100 M Q . A detailed discussion of 
the effects of varying the IMF on the properties of 
the present day elliptical galaxies will be presented 
in a forthcoming paper. 

The feedback due to supernova (SN) explosions 
moves the gas from the cold to the hot phase at a 
rate: 

M cold = --V(*)esN — , 

6 a 



with increasing M sp h, roughly as: 



(11) 



where tjsn is the number of Type II SNe expected 
per solar mass of formed stars (determined by 
the IMF, adopting a minimum progenitor mass of 
8 Mq), £?sn is the kinetic energy of the ejecta from 
each supernova (10 51 erg; e.g. Woosley & Weaver 
1986), and esN is the fraction of this energy which 
is used to reheat the cold gas. Analyses show that 
about 90% of the SN kinetic energy may be lost by 
radiative cooling (Thornton et al. 1998; Heckman 
ct al. 2000); we adopt esN = 0.05 as our reference 
value, consistent with the results by Mac Low & 
Ferrara (1999) and Wada & Venkatesan (2003). 

We relate V v [ r to the line-of-sight velocity dis- 
persion a using the relationship o~/V v - a ~ 0.65, 
derived by Ferrarese (2002) for a sample of local 
galaxies. 

The chemical evolution of the gas is followed by 
using classical equations and stellar nucleosynthe- 
sis prescriptions, as for instance reported in Ro- 
mano et al. (2002). 

We can elucidate the dependence of the star 
formation rate on halo mass by means of a sim- 
ple order of magnitude argument. Since the mean 
density within the virial radius is ~ 200/9 cr j t , the 
mean value of idyn is about a factor of 10 shorter 
than the expansion timescale, at all redshifts, in- 
dependently of the halo mass. If the gas is in virial 
equilibrium in the DM potential well, the effective 
cooling time of a pure hydrogen plasma is, assum- 
ing uniform density: 

i coo i !Off ~1.610 n (l + z vir r 5/2 /i 1/3 • 



(- 



\ 10 12 M, 



© 



1/3 



M vir /M g£ 
1/0.16 



1/2 



C" 1 yr(12) 



As discussed by Romano ct al. (2002, cf. their 
Fig. 10), for spheroidal galaxies with < 10 11 M 
in stars, the ratio between M vir and the mass in 
stars at z~0 (including remnants) M sph decreases 



M vir /M sph oc M t 



-1/3 
sph 



(13) 



due to the larger effect of stellar feedback in shal- 
lower potential wells. Since the mass of cooling gas 
is approximately equal to M sp h, i C ooi,cff turns out 
to be very weakly dependent on M v - lr . For an ef- 
fective clumping factor C a s > 10, Eq. (12), yields 
a cooling time shorter than the mean free-fall time 
for z > 3, for all galactic masses. Allowing for the 
gas metallicity obviously decreases i coo i- 

If the final mass in stars, -Ms pn , is proportional 
to the mass of the cooling gas, M spn , Eq. (13) 
implies that, for galaxies with M sp h < 10 11 M Q , 

M v i r cx oc Mgls- Thus the star forma- 

tion rate, ?/>(i) — M gas / max(i coo i, idyn) is approx- 

3 /2 

imately oc M vi ' r , since both timescales are ef- 
fectively independent of, or very weakly depen- 
dent on, M v j r . This means that stars form faster 
within larger dark-matter halos, as in the Anti- 
hierarchical Baryon Collapse scenario by Granato 
et al. (2001). If, as argued above, the star for- 
mation rate is controlled by idyn, we have ip(t) ~ 
32(M gas /10 n M Q )(l + z) 3 / 2 M yr" 1 . 

2.3. Black- hole growth 

As discussed by Haiman et al. (2003), the 
available information on the evolution of both the 
global star formation rate and the quasar emissiv- 
ity is broadly consistent with the hypothesis that 
star formation in spheroids and black-hole (BH) 
fuelling arc proportional to one another. 

One mechanism yielding such proportionality 
has been discussed by Umemura (2001), Kawakatu 
& Umemura (2002) and Kawakatu et al. (2003). 
In the central regions of proto-galaxies the drag 
due to stellar radiation may result in a loss of an- 
gular momentum of the gas at a rate that in a 
clumpy medium is well approximated by 



dlnJ 
dt 



L 



sph 



c 2 M v 



(1 



") 



(14) 



where £ spn is the global stellar luminosity and t 
is the effective optical depth of the spheroid. The 
latter quantity is given by r = T-/V; nt , where f is 
the average optical depth of single clouds and A^i nt 
is the average number of clouds intersected by a 
light ray over a typical galactic path. 
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The gas can then flow towards the center, feed- 
ing a mass reservoir around the BH at a rate 
(Kawakatu et al. 2003) 



inflow 



dlnJ 
3 dt 



J sph 



(i-O 



(15) 

During the early evolutionary stages the luminos- 
ity is dominated by massive main sequence stars, 
M > 5M Q , and is thus proportional to the star 
formation rate "4>(t). For the adopted IMF we 
have: 

M innow ~ 1.2xlO- 3 ^(t)(l-e- T ) Mgyr^ 1 . (16) 

While this expression is useful for the analytical 
estimates presented in Sect. 3.1, in the full calcu- 
lations we have adopted the values L spn computed 
using our spectrophotometric code GRASIL. We 
parameterize the optical depth r as: 



to 



Z 



1O 12 M 



(17) 



and for To we explore the range from 1 to 10 (see 
Sect. 3). 

An order of magnitude estimate of the relevant 
timcscales and luminosities can be derived assum- 
ing that -Mi n fi ow is directly the accretion rate on 
the central BH. Then, the timescale for it to grow 
to a mass Mbh is, using the order of magnitude 
estimate of il>(t) given in the previous section, 



^BH 



Mi 



BH 



M 



2.3 x 10 fc 



M, 



BH 



inflow 



10 8 M Q VH^Mt 



(1 + *) 



-3/2 



1 



-T\-l 

) yr ■ 



(18) 



Thus, if the accretion is not limited e.g. by ra- 
diation pressure or by angular momentum, black- 
holes can grow to very large masses from small 
(e.g., stellar mass) seeds in a time shorter than 
the age of the universe at all relevant redshifts. If 
Mbh oc M gas , the growth time is independent of 
M BH . 

For accretion with radiative efficiency 7/, the 
bolometric luminosity is L D oi — 77M nnow c 2 . Fol- 
lowing Elvis, Risaliti, & Zamorani (2002), we 
adopt, as a reference value, rj = 0.15. In units 
of the Eddington luminosity 

i Edd ~ 1.26 x 10 46 (M BH W/10 8 M Q )ergs- 1 , 

(19) 



and setting Mbh(^) — Mi nnow t, we have: 

-l 



Lbo\ _ , 3 rj 



L Edd 0.15 Vl0 8 yr 



t 



(20) 



showing that the fast accretion phase and the 
build up of super-massive BH can be quite fast, 
thus accounting for QSO at high z. 

Obviously, super-Eddington accretion is not the 
only mechanisms for a rapid black-hole growth. At 
the other extreme we may have accretion with low 
radiative efficiency such as the black-hole merg- 
ers at high redshifts advocated by Haiman et al. 
(2003), which may be testable by gravitational 
wave experiments like LISA (Mcnou, Haiman & 
Narayan 2001). Since these early evolutionary 
phases are expected to be heavily dust obscured, 
a powerful tool to discriminate among the various 
possibilities is hard X-ray emission. Far-IR/mm 
observations are also useful, but less capable of 
distinguishing the effect of the AGN from that of 
a starburst. 

While mechanisms for super-Eddington accre- 
tion have been proposed (Begelman 2001, 2002), 
fully unperturbed free-fall is unrealistic. Corre- 
spondingly, iboi/^Edd will not reach the extreme 
values indicated by Eq. (20) and the duration of 
the iboi/^Edd ^ 1 phase can be longer. In our 
model, we first let the material infalling at the 
rate given by Eq. (15) accumulate in a circum- 
nuclear mass reservoir, and then flow toward the 
black hole on a time scale depending on the vis- 
cous drag 

T v isc ~ I v . The reservoir accumulates 
mass at a net rate M res = M; nnow — Mbh, where 
Mbh is computed as follows. Following Duschl ct 
al. (2000) and Burkert & Silk (2001) we adopt a 
viscosity v = 1Z~ rit v r , where TZ C nt — 100-1000 is 
the critical Reynolds number for the onset of tur- 
bulence. With these assumptions the viscous time 
can be expressed as t v ; sc = Td yn 7?-crit.- The dynam- 
ical time Tdyn = (37r/32Gp s ) 1/2 is referred to the 
system 'black hole plus reservoir'. 

The accretion radius of the BH is given by 
r a = GMbh /Vyiy Defining the reservoir dimen- 
sion R Tes = ar ai we estimate p s as the mean 
density within a sphere of radius R res with mass 
M BH + M res , to get: 



Tdyn 



"(I) 



3/2 Q 



M 



3/2 
BH 



V v l (M BH + M rcs )V2- 



(21) 
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The viscous accretion rate onto the BH can be 
defined as 



MV1SC 
BH 



fcaccr G vM B J 



The constant fc„ 



M BH \ 

M res ; 

(22) 

7 r( a /2) 3 / 2 (Ki r /^) 3 ^crit]- 1 
has a rather wide range of possible values. For 
T^crit = 100 and a — 10 we have fc accr ~ 10" 4 that 
we adopt as the reference value in the following. 
The actual accretion rate is then given by 



1/2 



M BH = min (M^ C ,AM J 



Edd) 



(23) 



where M Edd = L Edd /(i]c 2 ) and A = (i/i E dd)max 
is the maximum allowed Eddington ratio. In our 
computations, we allow at most mildly super- 
Eddington accretion, i.e. A < a few. 

2.4. QSO feedback 

The QSO activity affects the interstellar medium 
of the host galaxy and also the surrounding in- 
tergalactic medium through both the radiative 
output and the injection of kinetic energy produc- 
ing powerful gas outflows. Quasar-driven outflows 
have been invoked to produce an intergalactic 
magnetic field (Furlanetto & Loeb 2001), and to 
preheat the intra-cluster medium (Valageas & Silk 
1999; Kravtsov & Yepes 2000; Wu et al. 2000; 
Bower et al. 2001; Cavaliere et al. 2002; Plata- 
nia et al. 2002). In the case of radio loud QSOs 
there is evidence that up to about half of the total 
power is in the jets (Rawlings & Saunders 1991; 
Celotti et al. 1997; Tavecchio et al. 2000). In 
broad absorption line (BAL) QSOs, the kinetic 
power of the outflowing gas can be a significant 
fraction of the bolomctric luminosity (Begelman 
2003). X-ray observations of BAL QSOs revealed 
significant absorption (Nh > I0 23 cm~ 2 ), imply- 
ing large outflows (M out > 5 M Q / yr) and large ki- 
netic luminosities Lk (Brandt et al. 2001; Brandt 
& Gallagher 2000). 

Theoretical studies on the mechanisms respon- 
sible for AGN-driven outflows show that efficient 
acceleration could be due to radiation pressure 
through scattering and absorption by dust (see e.g. 
Voit et al. 1993) and scattering in resonance lines 
(see e.g. Arav, Li & Begelman 1994). Murray et 
al. (1995) presented a dynamical model for a wind 
produced just over the disk by a combination of ra- 
diation and gas pressure. In a similar way Proga, 



Stone & Kallman (2000) showed that a wind can 
be launched from a disc around a 10 8 M black 
hole with velocities up to 0.1c and mass-loss rate 
of 0.5M©yr _1 . Following Murray et al. (1995), 
an approximate solution for the wind velocity pro- 
duced by line acceleration as a function of the ra- 
dius is: 

v = v oc (l- r J-) 2 ^ (24) 

where r/ is the radius at which the wind is 
launched. 



The asymptotic speed is 



^ GMB H y /2 



(25) 



where 7 is related to the force multiplier (see e.g. 
Laor & Brandt 2002). Adopting the reference val- 
ues of the parameters of Murray et al. (1995), we 
have 7 ~ 3.5. By replacing the BH mass with 
the corresponding Eddington luminosity £Edd,46 
in units of 10 46 ergs -1 , we get: 



6.2x10" 



10 



16, 



-1/2 



i E dd,4 6 1/2 . (26) 



Detection of outflows with velocities ranging be- 
tween 0.1c and 0.4c in the Broad Absorption Line 
(BAL) quasars APM08279 + 0522 (Chartas et al. 
2002) and PG1115 + 080 (Chartas et al. 2003) 
has been reported, based on observations of X-ray 
broad absorption lines performed with the Chan- 
dra and XMM-Newton X-ray observatories. In the 
case of APM08279 + 0522 an intrinsic bolometric 
luminosity of Lboi ~ 2.4 x 10 47 ergs _1 has been 
estimated by Egami et al. (2000). Under the as- 
sumption of L E dd ~ 1"3 Z/boh the maximum ob- 
served velocity v ~ 0.4c is obtained with a launch- 
ing radius 77 ~ 0.5-2 x I0 16 cm. Such small val- 
ues of rj are confirmed by the observed variabil- 
ity of the absorption line energies and widths in 
APM08279 + 0522 over a proper timescale of 1.8 
weeks (Chartas et al. 2003). 

The asymptotic speed is reached at r > 40ry. 
If f c is the covering factor of the outflow and using 
M w = 4irr 2 p(r)v r - 47r/ c m ff iVH40r / i; oo , we get 



M w = 2.6/^4^ 



1/2 



10 16 Cl!: 

Ah/10 22 



where and N 2 2 — iV#/lU"cm 
1 .5x 10 16 cm as a reference value, the kinetic power 



M yr 1 
(27) 

; Adopting rj — 
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in the outflow is 



removes it from the cold phase at a rate 



L K = -M w vl ~ 3.6 X 10 44 /c^224dd46 er g S " 



M. 



QSO 



2 L h M, 



cold 



(28) 

and may thus amount to several percent of the ac- 
cretion luminosity for highly luminous QSOs. In- 
terestingly, this corresponds to the power required 
to account for the pre-heating of the intra-cluster 
medium (Bower et al. 2001; Platania et al. 2002; 
Lapi et al. 2003). It is also interesting to notice 
that M^/Macc 2.5(e/0.1)/ c 7V 22 i E ^ 2 46 , imply- 
ing that for highly luminous QSOs M w ~ M acc , 
if f c N 2 2 ~ 1. High luminosity QSOs emitting at 
the Eddington limit are able to generate winds in- 
volving relatively small amounts of gas, but with 
very high velocities and significant kinetic ener- 
gies. When the luminosity decreases below the 
Eddington limit we replace the Eddington lumi- 
nosity with the bolometric luminosity in Eqs. (28) 
and (30). 

Estimating the fraction of the kinetic luminos- 
ity Lk transferred to the interstellar medium is 
a rather complex problem. One possible effect of 
outflows and jets is to transport the ambient gas 
to larger radii. This effect has been recently esti- 
mated through numerical simulations by Bruggcn 
et al. (2002). By investigating the effects of radio 
cocoons on the intra-cluster medium, they con- 
cluded that frequent low-power activity cycles are 
rather efficient in stirring the environment gas, 
particularly in regions close to the injection point 
of the high speed gas. Enfilin & Kaiser (2000) 
evaluated the energy accumulated in the cocoons 
around radio-galaxies. Inoue & Sasaki (2001) ar- 
gued that for non relativistic cocoon plasma, the 
fraction of the jet energy deposited in the intra- 
cluster medium can reach 40%, neglecting the ra- 
diative cooling. Similar conclusions have been 
reached by Bickncll et al. (1997). Nath & Roy- 
chowdhury (2002), taking into account also the ra- 
diative losses, found that a large fraction (> 0.5) 
of the kinetic energy of BAL and radio loud QSOs 
is transferred to an ambient gas with number den- 
sity n ~ 0.1-1 cm -3 and temperature T v i r ~ 10 6 
K. These values of density and temperature are 
quite similar to those of the gas in the outer re- 
gions of massive galactic halos. 

We assume that the QSO feedback heats up 
the interstellar medium at a rate Lh = fhLx and 



cold 



3 cr 2 M s 



(29) 



where M gas = M cold 



Mj n f is the mass of the gas 
in the cold phase plus that of the gas which has 
not yet fallen in the star forming region. Setting 
eqso = (/fc/0.5)(/ c /0.1)(JV 22 /10), we get 



r 3/2 

ON/in 3 £ QSO^Edd,46 Mcold _i 

- 2Xl ° ( g /300kms-i) 2 M gas Meyr ' 

(30) 

We have explored the range 1 < £qso < 10, adopt- 
ing eQso = 6 as our reference value. The fraction 
-Mcold /M gas in the right-hand side of Eq. (29) has 
been introduced in order to share the QSO feed- 
back between the cold and the infalling gas. Thus 
the QSO feedback also removes the infalling gas 
at a rate 

,v qso _ _ 2 Lh M inf 



M, 



inf 



3 (T 2 M* 



(31) 



The quasar feedback can easily heat the inter- 
stellar gas to temperatures ~ IkeV (Valageas 
& Silk 1999; Bower et al. 2001; Nath & Roy- 
chowdury 2002; Platania et al. 2002; Lapi et 
al. 2003), thus unbinding it and making it flow 
into the intcrgalactic medium (IGM). The cor- 
responding energy injection in the IGM may ac- 
count for the steepening of the X-ray luminosity- 
temperature correlation observed in groups of 
galaxies (O 'Sullivan, Ponman, & Collins 2003). 
Given the low IGM gas density, only a small frac- 
tion of it will eventually cool down and fall back 
again to form stars. 

Table 1: Model parameters, reference value and 
number of the equation where it is introduced (or 
short description). 
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Fig. 1. — Scheme of the baryonic components in- 
cluded in the model (boxes), and of the corre- 
sponding mass transfer processes (arrows). The 
numbers near the arrows point to the main equa- 
tions describing those processes. 

3. Results 

The results presented here refer to halos with 
mass in the range 2.5 x 10 11 M Q < M vir < 1.6 x 
10 13 Mq, formed at z v ; r > 1.5. This lower limit to 
the virialization redshift allows us to crudely fil- 
ter out halos hosting discs and irregular galaxies. 
In our view, late-virialized objects had more time 
to acquire a substantial angular momentum and 
are less likely to take on an early-type morphol- 
ogy. This is in keeping with the common notion 
that stellar populations in disks are on average 
substantially younger than those in spheroids. 

The model, whose main baryonic components 
and corresponding mass transfer are sketched in 
Fig. 1, allows us to follow the star formation and 
chemical enrichment histories of spheroidal galax- 
ies, as well as the growth of their central BHs, 
once the virialization time and the halo mass is 
given. Interfacing it with the code GRASIL (see 
Silva et al. 1998 for details), which computes the 
spectrophotometric evolution from radio to X-ray 
wavelengths including the dust effects, we get the 
spectral properties of the spheroidal galaxies as 
a function of cosmic time. In the application of 
GRASIL, we have made the standard assumption 
that the starburst is highly obscured throughout 



its duration, with an 1 /im optical depth of molec- 
ular clouds, where new stars are born, of n = 30 
for solar metallicity, and scaling linearly with Z. 
As for the dust emissivity index, we have adopted 
the canonical value of 2. The results presented in 
this paper are only weakly sensitive to variations 
of n by a factor of several. 

The relationship between n and the optical 
depth given by Eq. (17) is not straightforward. 
The starlight is processed by dust within molec- 
ular clouds and re-emitted in the far-IR, with a 
peak at 60-100 ^m. At these wavelengths, the op- 
tical depth is lower than that at 1 /xm by a factor 
~ 10~ 2 . The effective optical depth to the inter- 
cloud radiation, responsible for the drag exerted 
on gas clouds as discussed in Sect. 2.3, is there- 
fore r ~ 10 _2 TiiVi n t, where N[ nt is the mean num- 
ber of clouds intersected by a photon (Kawakatu 
& Umemura 2002). The latter number depends, 
among other things, on the radial distribution of 
the clouds within the galaxy. For instance, a pop- 
ulation of giant molecular clouds with a typical 
size of 10 pc, masses of the order of 10 5 M Q , uni- 
formly distributed in the central ~ 10 kpc of the 
galaxy, give N int ~ 10, for M gas ~ 10 12 M Q . As 
anticipated, for t [Eq. (17)] we explored the range 
1-10. 

In this Section we present some of the most im- 
portant results, obtained using the reference val- 
ues of the parameters already discussed and sum- 
marized in Table 1 . The discussion of the effects of 
varying the parameters within the allowed ranges 
will be presented in the next Section. 

Examples of the time evolution of the various 
components for different virialization redshifts and 
halo masses are shown in Fig. 2 and Fig. 3. 

3.1. Time delay between star formation 
and QSO activity 

A key result of the model is the prediction of 
the time delay between the onset of vigorous star 
formation, at the virialization epoch, and the peak 
of the QSO activity. The combined action of stel- 
lar and nuclear feedbacks, eventually sweeping out 
the interstellar medium, determines the duration 
of the star formation. 

As shown in Fig. 4 the duration of the most ac- 
tive star formation phase decreases with increasing 
M v i r and z v - lr . The timescale Tioo is defined as the 
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Fig. 2. — Evolution of the mass in stars, in the cold 
(M co id) and infalling (Mj n f) gas components, and 
of the Star Formation Rate, as a function of the 
galactic age T, for 3 values of M v - lr and z v i r = 4. 
Note the sharp break corresponding to the sweep- 
ing out of the interstellar medium by the quasar 
feedback, for large M vir . 



galactic age at which the gas mass still available 
for star formation or accretion, i.e. M coo i + M;„f , 
is reduced to 1% of the initial value. This corre- 
sponds approximately to the duration of the star 
formation burst empirically estimated by Granato 
et al. (2001), and given by their Eq. (8). In the 
redshift range 3 < z y -„ < 6, Ti o — 0.5-1 Gyr 
for M v ; r > 10 12 M Q , corresponding (see Fig. 5) to 
M* P h > 2.5 x 10 10 Mq. For smaller masses the du- 
ration of the actively star forming phase increases 
to reach 1.5-3 Gyr for M sph ~ 3 x 10 9 M , if virial- 
ized at z v -„ > 3, and is even longer if virialization 
occurs at lower z. 

A feeling of the role of the main ingredients of 
the model in shaping this result can be obtained 
ignoring the reservoir, so that the inflow onto the 
central black hole is governed by Eq. (15). The 
inflow is definitively halted at a time t out when 
the energy injected by the QSO in the interstellar 



Fig. 3. — Growth of the black-hole and of the 
reservoir as a function of the galactic age T, for 
the same values of M v i r and z v ; r = 4 as in Fig. 2. 
The apex of dMBn/dt corresponds to the end of 
the exponential growth of the black- hole mass, 
when the mass accretion rate becomes insufficient 
to support the adopted maximum Eddington ra- 
tio, (L/IzEdd)max = ^- 

medium equals the binding energy of the gas (cf. 
Cavaliere et al. 2002): 



tout 



f h L K dt = M vil V? 



2 

vir ' 



(32) 



where VW is related to the mass inside the virial 
radius, M vil} by (Bullock et al. 2001b): 



1/2 



M ■ \ 1/3 

Jvi vir \ i _1 

kms 



Kir = 75(1 + Z vir ) 

(33) 

Using Eq. (28) and Eq. (18) with cxp(-r) < 1, 
and assuming that the black hole growth occurs 
in a time shorter than the expansion timescale, 
we obtain: 



i out ~ 7.9xl0 8 h i / 15 f- 2/5 (l + z vil )- 1 
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Fig. 4. — Duration of the star-formation burst as 
a function of M V - 1T , for 3 values of z v ; r (see text 
for the exact definition of Tioo)- For our choice of 
the parameters, at low z v i r the cooling time be- 
comes longer than the dynamical time at large 
masses, causing the upturn of Tioo- The long- 
dashed line close to the solid line corresponding to 
z vir = 4 represents the empirical recipe adopted 
by Granato et al. (2001). 

~ 1.2 x 10 9 ft 4 / 15 /fc 2/5 ■ 

/ M ■ \ ~ 7/3 ° 
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where the last step follows from Eq. (13), normal- 
ized to M vir /M sph ~ 20 for M sph = 5 x 10 10 M Q , 
assuming M gas ~ M sph . 

The final black-hole mass (except for the effect 
of later re-activation phases) can be roughly esti- 
mated to be ~ AfBH(^out): 

M BH (Ut) * 4.3 x 10W 16 K 2/S ■ 
( _M^_\ 2/3 ( M gas V /5 

~ 3.2 x 10 7 h^ 15 f- 2/5 ■ 
( M ■ \ 19/15 

/ \ 19/5 /-. \ -9/10 



Fig. 5. — Virial to stellar mass ratio at the present 
time as a function of M v i r for 3 values of z v -„. 

where we have used Eqs. (13) and (33), and have 
adopted the relationship a ~ 0.65Kir (Ferrarese 
2002). This result is in remarkably good agree- 
ment with the determination by Gebhardt et al. 
(2000) of the Mbh-c relationship. A full compar- 
ison of our model predictions with observational 
data is shown in Fig. 6, where the expected spread 
due to the distribution of virialization rcdshifts is 
also illustrated. Note the predicted fall-off of Mbh 
at low values of a, due to the combined effect of 
SN feedback which is increasingly efficient with de- 
creasing halo mass in slowing down the gas infall 
onto the central BH, and of the decreased radia- 
tion drag due to a decrease of r [Eq. (17)], which 
is, for these objects, -C 1. This steepening of the 
lower part of the Mbh-c relationship translates 
in a flattening of the lower part of the Mv^-Mbh 
correlation, in agreement with the results reported 
in Fig. 5 of Ferrarese (2002). 

In Fig. 7 the predicted BH mass function at 
z = is shown against the local BH mass func- 
tion recently derived by Shankar et al. (in prepa- 
ration) following Salucci et al. (1999). The 
total local BH mass density is estimated to be 
5.2 x lO 5 M Mpc~ 3 . 

The physics that rules the fast growth of the 
central BH in massive halos is the key point: its 
growth is paralleled by an increasing feedback 
from the nuclear activity. In the most massive ha- 
los the quasar feedback eventually removes most 
of the gas (and dust), leaving the nucleus shining 
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Fig. 6. — Predicted relationship between black- 
hole mass and line-of-sight velocity dispersion of 
the host galaxy for different virialization redshifts. 
Data are from Tremaine et al. (2002). 

as an optical quasar until the reservoir mass is ex- 
hausted on a timescale ~ 10 7 yr (see Fig. 3). Since 
then the host galaxy evolves passively and the BH 
becomes dormant (apart from possible reactiva- 
tions). The process slows down with decreasing 
halo mass, as the supernova feedback becomes in- 
creasingly important. 

3.2. Photometric properties and metal 
abundances of spheroidal galaxies 

In small halos (M vir < 10 11 M ) the star for- 
mation rate is kept low by stellar feedback, which 
heats up most of the gas and moves it to lower den- 
sity outskirts, where the cooling time is very long. 
Only when the virial mass exceeds a few 10 11 M 
the potential wells are deep enough to allow a 
more effective star formation. Figure 5 shows the 
dependence on M v ; r of the ratio Mvi r /M sprl , for 
three values of z v ; r . Bright galaxies virializing at 
z < 4 are predicted to have M v i r /M spri — 40. 
For comparison, McKay et al. (2002) report, for 
bright galaxies, M/L values in the SDSS (/-band of 
171 ± 40 based on the dynamics of satellites, and 
of 270 ± 35 based on weak lensing measurements; 
adopting for the stellar component of early-type 



galaxies M/L v ~ 4-5.9 (Fukugita et al. 1998) this 
translates (for Ly ~ L g ) into M v i r /M sp h — 30-70. 

Once the star formation law is specified, the 
chemical evolution of the galaxy is followed by us- 
ing classical equations and stellar nucleosynthe- 
sis prescriptions (Granato et al. 2001). The pre- 
dicted mean metallicity, (Z+) in stars is shown, as 
a function of the mass in stars, in Fig. 8. The 
average metallicity increases from about half so- 
lar to super-solar with increasing galaxy mass. In 
particular, for our reference model, the metal con- 
tent is super-solar for M v - K > 5 x 1O 11 M , or 
M sph >10 10 M Q . 

It should be stressed that this result refers to 
the average metallicity of the galaxies, since our 
model is one-zone. In the central regions, star- 
formation, and consequently metal enrichment, 
is faster because of the higher densities entail- 
ing shorter dynamical times. Adopting densities 
appropriate for the innermost 10% of the galaxy 
mass, the predicted metallicity of the gas is be- 
tween 3 to 5 times solar, with a general trend at 
increasing with M v ; r . This result compares fairly 
well with several recent direct chemical measure- 
ments of circum-nuclear gas in high-z QSO (e.g. 
Dietrich et al. 2003). 

Another evident trend is the increase with the 
galaxy mass of the ratio between the abundance 
of the so-called a-elements and that of Fe. Fig. 9 
shows the mean [Mg/Fe] ratio, in stars, as a func- 
tion of the galaxy mass for three different values 
of the virialization rcdshift. The ratio [Mg/Fe] 
is always super-solar, and is further enhanced 
for M v j r <; 10 12 M Q . However, in spite of the 
high metallicity reached by our models, we find 
that the average stellar [Fe/H] is always subso- 
lar (< —0.25). This trend is dictated by the mass 
dependence of Tioo: in high mass objects the star- 
formation activity is essentially ended when SNIa 
start to pollute the ISM with Fe. 

Both the trend of the global metallicity and the 
ratio of the a- elements to Fe abundance with 
the galaxy mass, are in fair agreement with what 
is inferred from the spectrophotometric observa- 
tions of local early type galaxies. The narrowness 
and the inclination of the color-magnitude rela- 
tion of elliptical galaxies hints at a rapid formation 
process and at a range in metallicity from about 
half solar to twice solar (e.g. Bressan, Chiosi & 
Fagotto 1994). Furthermore, narrow band obser- 
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Fig. 7. — Predicted local black-hole mass function 
(solid line) compared with the recent estimate by 
Shankar et al. (in preparation, region delimited by 
the two thin dashed lines) . The total mass density 
in BHs they derive is p B H - 4 x lO 5 M Mpc~ 3 , 
25% less than our model. The decline of the model 
at low Mbh is due to having considered only ob- 
jects with M vir > 2.5 x lO n M0. 



vations not only provide evidence for a similar 
spread in metallicity (Bernardi et al. 2003b), but 
also demand a significant enhancement of the a- 
elements in the more massive ellipticals (Worthcy 
et al. 1994). The observed run of Mg and Fe 
narrow band indices with the luminosity of the 
galaxy or with the H(3 index, cannot be repro- 
duced by stellar population models based solely 
on solar partition of heavy elements (Worthey et 
al. 1994; Trager et al. 2000a,b; Bernardi et al. 
2003b). 

It is worth noticing at this level that while the 
model reproduces, for a relatively broad range of 
parameter choices, the main trends of the metal 
abundances inferred from observations of local el- 
lipticals, our predictions cannot yet be directly 
compared with observations. For example, nar- 
row band indices show a complicated dependence 
on global metallicity, partition of heavy elements 
and age of the stellar populations, that render the 
adoption of common scaling relations (e.g. Mat- 
teucci, Ponzone & Gibson 1998) quite uncertain. 
We should thus generate mock catalogues of lo- 
cal early type galaxies and investigate their spec- 
trophotomctric properties by means of adequate 
spectral synthesis tools. Work in this direction is 
in progress and will be reported in a subsequent 
paper. 
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Fig. 8. — Mean stellar metallicity as a function of 
the present-day stellar mass for the usual 3 values 
of z vir . 
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In these galaxies the K-band luminosity is a 
quite good indicator of the mass in stars. In 
Fig. 10 we compare the modelled K-band local 
luminosity function with the recent determina- 
tions by Kochanek et al. (2001) and Huang et 
al. (2003). As already remarked, the solid line in- 
cludes the contribution of all spheroids with mass 
in the range 2.5 x 10 11 < M vir < 1.6 x 10 13 , formed 
at Zvir ^ 1-5. According to our model, spheroidal 
galaxies born in halos with M vir < 2.5 x 10 11 M Q 
have masses in stars M sp h < 10 9 M© and K-band 
magnitudes > —21.3 (cfr. Fig. 5). 

For a more direct comparison with the observa- 
tional data, the dot-dashed line in Fig. 10 shows 
the predicted local luminosity function for early- 
type galaxies only. It was obtained subtracting 
from the results shown by the solid line, the contri- 
butions of bulges of Sa and Sb galaxies estimated 
using the local B-band luminosity functions and 
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the bulge-to-total luminosity ratios for galaxies of 
these morphological types derived by Salucci et al. 
(1999) and adopting a color B-K = 4.1. It should 
also be noted that the mass in stars for small ha- 
los is strongly dependent on the stellar feedback, 
namely on the IMF and on the SN efficiency csn: 
an increase of the efficiency from 0.1 to 0.3 results 
in a decrease of mass in stars by a factor ~ 3. This 
dependence weakens with increasing M v ; r . 

Figure 11 compares the predicted rest-frame 
K-band luminosity function at z = 1.5 with 
the observational determination by Pozzetti et 
al. (2003), which, at the highest luminosities, is 
substantially above predictions of the models by 
Kauffmann et al. (1999), Cole et al. (2000), and 
Menci et al. (2002). 

3.4. SCUBA galaxies and EROs 

In large galactic halos the star formation rate 
turns out to be very high at high redshift, yield- 
ing a quick increase of the metallicity and of the 
dust mass. The latter is computed by GRASIL 
as proportional to the product of the gas mass 
by its metallicity, with a coefficient determined by 
the condition of a gas-to-dust ratio of 110 for so- 
lar metallicity (Silva et al. 1998). Thus, most of 
the star formation occurs in a dusty environment, 
so that these galaxies are powerful far-IR/sub- 
mm sources, highly obscured in the visual and 
near IR bands. In Fig. 12 we have plotted the 
model predictions at 850 (im against the SCUBA 
counts. The predicted redshift distributions are 
similar to those shown by Granato et al. (2001). 
In particular, for a flux density limit of 5 mJy, 
the model gives a median redshift of 2.2 and an 
inter-quartile range of 1.6-3.3, to be compared 
with ^median = 2.4 and the inter-quartile range 
of 1.9-2.8 found for the sample of Chapman et al. 
(2003). 

After the interstellar medium has been swept 
out, galaxies evolve passively. The combination of 
redshift and aging soon makes them extremely red. 
We computed the expected contribution to the ex- 
tragalactic K-band counts of spheroidal galaxies 
in this phase. A comparison (Fig. 13) of the pre- 
dicted with the observed redshift distribution of 
galaxies with K < 20 (Cimatti et al. 2002b) shows 
that they fully saturate the high redshift tail of the 
distribution. 



Cimatti et al. (2002a) selected a complete sam- 
ple of Extremely Red Objects (EROs) {(R-K a ) > 
5) with K < 19.2. More than 60% of the ob- 
jects have redshift, mostly spectroscopic. On the 
basis of the spectra, the sample has been subdi- 
vided into dusty and non-dusty EROs. Their data 
suggest that there is a significant number of old 
dust-free ellipticals in place at z> 1. This result is 
confirmed by the subsequent analysis by Pozzetti 
et al. (2003), who found that the bright end of 
the K-band luminosity function at z > 1 is dom- 
inated by red/early type galaxies. Our model is 
consistent with these results (cf. Fig. 11). 

Our expectation of the existence of a significant 
population of luminous red galaxies at substantial 
redshifts is also borne out by the results of the very 
deep near-infrared photometry of the Hubble Deep 
Field-South with ISAAC on the VLT (Franx et 
al. 2003) and by follow-up Keck spectroscopy (van 
Dokkum et al. 2003). Model predictions closely 
match the surface density (3 ± 0.8arcmin~ 2 ) and 
the median redshift (2.6) of galaxies with K s < 
22.5 and ( J s - K s ) > 2.3 (Franx et al. 2003). 

4. Discussion 

To highlight the dependence of galaxy proper- 
ties on their halo masses, we make reference to 
three characteristic values, namely M v - lr ~ 10 114 , 
10 12 - 4 and 10 13 - 4 , referred to as low, intermediate 
and high masses, respectively. 

The clumping factor C affects only the evolu- 
tion of high and intermediate masses, where a de- 
crease of C from 20 to 1 strongly inhibits the ca- 
pability of gas to quickly cool down, form stars 
and feed the BH. For large galactic masses, the 
star-formation activity and the growth of the BH, 
rather than being confined to the first ~ 1 Gyr 
after virialization, continues for a time compara- 
ble to the Hubble time. However, the final ratio 
A^BH/Mjph is almost the same, both masses ul- 
timately depending on the SFR. Conversely, at 
lower masses the collapse becomes increasingly 
limited by the dynamical time even when C = 1. 

The supernova efficiency £sn has some effect on 
the evolution at all masses: M sp h, < > and 
MBn/Msph decrease with increasing esn, while 
Tioo increases. Differences are modest at high 
masses, but become dramatic at low masses. At 
M vir < 2 x 10 11 M0 and e SN > 0.2, the evolu- 
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tion of the star-formation rate can show damped 
oscillations at early epochs (T < 1 Gyr) . 

By converse, the QSO efficiency, cqso, influ- 
ences only intermediate and large masses, where 
an increase of this parameter determines a shorter 
active star forming phase, and lower values of (Z*), 
M sp h and Mqso, while the ratio Mbh/MspIi is lit- 
tle affected. An increase of the radiative efficiency 
e has a very similar effect. 

When the accretion is Eddington-limited (or, 
at most, mildly super-Eddington) , Tioo, {%*)■> and 
decrease with increasing seed BH mass: less 
e-folding times are necessary to produce substan- 
tial effects on the environment. Typically, M sp h 
decreases by 30% when the seed mass is increased 
from 10 3 to 10 4 M Q , but the ratio M B H/M sph de- 
creases by less than 10%. As already remarked, 
the development of low masses is instead weakly 
affected by the BH feedback. By converse, if sub- 
stantially super-Eddington accretion is allowed, 
the model is quite insensitive to the precise choice 
of the seed BH mass. 

The ratio between QSO and SN feedback is 
an increasing function of the virial mass of the 
galaxy (see Fig. 14). For low mass galaxies the 
integrated effect of the QSO is almost negligible 
compared to that of SNae, but takes over (typi- 
cally by a factor of a few) for intermediate masses 
(M v i r ~ 10 12 4 M©), and dominates (by a factor 
> 10) at high masses (M vir - 10 13 - 4 M Q ). Note 
that the QSO effect usually increases exponen- 
tially with time, while that of SNae increases more 
slowly. Thus the instantaneous QSO effect be- 
comes dominant, if ever, only a few e-folding times 
before the maximum of QSO activity. 

As mentioned at the beginning of Sect. 3, our 
definition of a galactic halo associated with a 
spheroidal galaxy is rather crude. However, the 
successful comparison with the observed popu- 
lation properties of spheroidal galaxies (Figs. 7, 
10, 11, 12, 13, and 15) gives us some confidence 
about the meaningfulness of the criterion we have 
adopted. 

Although we have not addressed the details of 
the formation of disk (and irregular) galaxies, we 
envisage them as associated primarily to halos viri- 
alizing at z v i r < 1.5, which have incorporated, 
through merging processes, a large fraction of ha- 
los less massive than 2.5 x 10 11 Mq virializing at 



earlier times. These low mass halos virialized at 
early times may become the bulges of late type 
galaxies. 

5. Summary and conclusions 

We have presented a detailed, physically grounded, 
model for the early co-evolution of spheroidal 
galaxies and of active nuclei at their centers. The 
model is based on very simple recipes, that can be 
easily implemented. In summary, we start from 
the diffuse gas within the dark matter halo falling 
down into the star forming regions at a rate ruled 
by the dynamic and the cooling times. Part of 
this gas condenses into stars, at a rate again con- 
trolled by the local dynamic and cooling times. 
But the gas also feels the feedback from super- 
novae and from active nuclei, heating it and pos- 
sibly expelling it from the potential well. Also the 
radiation drag on the cold gas decreases its angu- 
lar momentum, causing an inflow into a reservoir 
around the central black hole. Viscous drag then 
causes the gas to flow from the reservoir into the 
black hole, increasing its mass and powering the 
nuclear activity. 

In the shallower potential wells (corresponding 
to lower halo masses and, for given mass, to lower 
virialization redshifts), the supernova heating is 
increasingly effective in slowing down the star for- 
mation and in driving gas outflows, resulting in an 
increase of star/dark- matter ratio with increasing 
halo mass. As a consequence, the star formation 
is faster within the most massive halos, and the 
more so if they virialize at substantial redshifts. 
Thus, in keeping with the proposition by Granato 
ct al. (2001), physical processes acting on baryons 
effectively reverse the order of formation of galax- 
ies compared to that of dark-matter halos. 

A higher star-formation rate also implies a 
higher radiation drag, resulting in a faster loss of 
angular momentum of the gas (Umcmura 2001; 
Kawakatu & Umemura 2002; Kawakatu et al. 
2003) and, consequently, in a faster inflow towards 
the central black-hole. In turn, the kinetic energy 
carried by outflows driven by active nuclei through 
line acceleration is proportional to M^Q (Murray 
et al. 1995), and this mechanism can inject in the 
interstellar medium a sufficient amount of energy 
to unbind it. The time required to sweep out the 
interstellar medium, thus halting both the star for- 
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mation and the black-hole growth, is again shorter 
for larger halos. For the most massive galaxies 
(M vir > 10 12 M Q ) virializing at 3 < z v ir < 6, 
this time is < 1 Gyr, so that the bulk of the star- 
formation may be completed before type la super- 
novae can significantly increase the Fe abundance 
of the interstellar medium; this process can then 
account for the a-enhancement seen in the largest 
galaxies. 

The interplay between star formation and nu- 
clear activity determines the relationship between 
the black-hole mass and the mass, or velocity dis- 
persion, of the host galaxy, as well as the black- 
hole mass function. As illustrated by Figs. 6 and 
7, the model predictions arc in excellent agree- 
ment with the observational data. A specific pre- 
diction of the model is a substantial steepening 
of the Mbh~ct relation for a < 150 kms -1 : the 
mass of the BH associated to less massive halos is 
lower than expected from an extrapolation from 
higher masses, because of the combined effect of 
SN heating, which is increasingly effective with de- 
creasing galaxy mass in hindering the gas inflow 
towards the central BH, and of the decreased ra- 
diation drag [see Eq. (15), with t<1]. 

Coupling the model with GRASIL (Silva et al. 
1998), the code computing in a self-consistent way 
the chemical and spectrophotometric evolution of 
galaxies over a very wide wavelength interval, we 
have obtained predictions for the sub-mm counts 
and the corresponding redshift distributions as 
well as for the redshift distributions of sources de- 
tected by deep K-band surveys, which proved to 
be extremely challenging for all the current semi- 
analytic models. The results, shown by Figs. 12 
and 13, are again very encouraging. 

A discussion of the evolutionary properties of 
AGNs predicted by the present model is deferred 
to a future paper, where we will address, among 
other things, the complex issue of how the bolo- 
metric luminosity produced by accretion processes 
shows up in different electromagnetic bands. We 
note however that the analysis by Granato et al. 
(2001), who approached the problem the other 
way round, i.e. inferred the formation history of 
spheroidal galaxies and of galactic bulges from the 
observed epoch-dependent luminosity function of 
quasars, got results in detailed agreement with 
those presented here. Thus an at least qualita- 
tive agreement of the present model with the data 



on AGNs seems to be ensured. This is confirmed 
by the successful comparison of our model pre- 
dictions for the relationship among the black-hole 
mass and the velocity dispersion of the host galaxy 
(Fig. 6) , and for the local black- hole mass function 
(Fig. 7). Also, the predicted history of global ac- 
cretion rate onto the central black-holes (Fig. 15), 
which, in our scheme, is directly proportional to 
the history of the bolometric luminosity density 
produced by AGNs, with its peak in the redshift 
range 2 to 3, is nicely consistent with the results 
of optical surveys (see, e.g., Fig. 8 of Fan et al. 
2003). 
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Fig. 9. — Mean Mg/Fe abundance ratio as a func- 
tion of the present-day stellar mass for the usual 
3 values of z v ; r . 




Fig. 10. — Predicted K-band local luminos- 
ity function of massive spheroids (solid line) 
compared with observational determinations by 
Kochanek et al. (2001; diamonds) and Huang et 
al. (2003; stars). The dot-dashed line shows the 
predicted local luminosity function of early-type 
galaxies only, i.e. after subtracting from the global 
luminosity function of spheroids the contributions 
of bulges of Sa and Sb galaxies (see text). 

This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.0. 
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Fig. 11. — Predicted K-band luminosity func- 
tion of massive spheroids at z — 1.5 compared 
with observational determination by Pozzetti et 
al. (2003). 
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Fig. 13. — Predicted redshift distribution of galax- 
ies brighter than K = 20 compared with the re- 
sults of the K20 survey (Cimatti et al. 2002b). 
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Fig. 12.— Predicted 850 {im extragalactic counts 7 7 5 8 8 5 9 9 5 

compared with SCUBA counts by Blain et al. ° 9 y ' 

(1999, open circles), Hughes et al. (1998, star 

and triangles), Barger et al. (1999, open squares), Fi S- 14 -— Ratio between the mass ejected from 

Eales et al. (2000, filled circles), Chapman et al. the galaxy per unit time by the feedback of the 

(2002, filled triangle), and Borys et al. (2002, filled Q S0 and by that of SNae. 
squares). 
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Fig. 15. — Predicted star formation in spheroids 
and BH mass accretion rates per unit volume as a 
function of rcdshift. 
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